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Abstract:  The  Tomahawk  land-attack  cruise  missile  can  be  launched  from  a  ship  or  submarine,  and  can  deliver  its  warhead 
precisely  to  a  target  at  long  range.  There  are  several  variants  of  the  Tomahawk  missile,  each  with  specialized  capabilities.  For  each 
Tomahawk  Missile  Sequence  Number  (MSN)  task  (i.e.,  mission),  the  Tomahawk  missile  variants  can  be  ranked  with  respect  to 
their  ability  and  cost  effectiveness  to  perform  that  task.  A  given  land  attack  strike  order  may  involve  a  large  number  of  Tomahawk 
missiles  and  numerous  Tomahawk  launch  platforms.  Operational  planners  select,  in  real  time,  feasible  launch  platforms  to  execute 
Tomahawk  taskings.  The  Tomahawk  tasking  in  a  strike  order  includes  not  only  primary  assignments  but  also  backup  assignments.  On 
board  each  launch  platform,  the  precise  allocation  of  specific  Tomahawk  missiles  to  the  Tomahawk  MSN  task  assigned  is  optimized 
with  a  model  described  here.  We  help  naval  operational  planners  select,  in  real  time,  appropriate  feasible  launch  platforms  to  fulfill 
the  Tomahawk  tasking  in  a  strike  order.  ©  2009  Wiley  Periodicals,  Inc.  Naval  Research  Logistics  58:  281-295,  201 1 
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“Be  swift  like  lightning  in  the  execution; 

And  let  thy  blows,  doubly  redoubled, 

Fall  like  amazing  thunder  on  the  casque 
Of  thy  adverse  pernicious  enemy” 

Richard  //,  Shakespeare 

1.  INTRODUCTION 

The  Tomahawk  (Fig.  1)  is  a  long-range,  subsonic,  land- 
attack  missile  manufactured  by  Raytheon.  It  is  a  tactical, 
conventional  weapon  that  can  be  launched  from  U.S.  Navy 
surface  ships  and  submarines,  and  some  Royal  Navy  sub¬ 
marines.  It  flies  a  planned  route,  delivering  its  warhead  pre¬ 
cisely  to  a  target  at  long  range,  at  a  specified  time.  A  given 
land  attack  strike  order  may  involve  a  large  number  of  Toma¬ 
hawk  missiles  of  different  types  (called  Tomahawk  variants) 
and  numerous  Tomahawk  launch  platforms. 

In  this  article,  we  describe  a  “predesignation”  optimization 
model  and  a  fast  heuristic  algorithm  to  help  the  U.S.  Navy’s 
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operational  planners  select,  in  real  time,  feasible  launch  plat¬ 
forms  to  fulfill  the  tasking  in  a  strike  order.  Predesignation  is 
the  process  by  which  a  set  of  one  or  more  specific  Tomahawk 
missiles,  on  one  or  more  specific  firing  platforms  is  matched 
to  the  mission  requirements  in  a  “strike  plan.”  For  a  given 
target  (that  is,  a  set  of  aim  points  to  be  attacked  in  a  single, 
coordinated  action),  a  strike  plan  specifies  how  many  and 
what  kind  of  Tomahawk  missiles  are  requested.  In  addition, 
we  address  specific  aspects  of  a  goal  hierarchy  involved  in 
the  coordination  of  a  large  number  of  Tomahawk  missiles 
and  naval  platforms. 

There  is  a  rich  and  long-standing  literature  on  missile 
assignment.  For  example,  the  early  work  of  Bracken  and 
McCormick  [4]  presents  a  nonlinear  approximation  to  a 
weapon  allocation  problem  where  the  goal  is  to  maximize 
expected  inflicted  damage.  Matlin  [26]  reviews  much  of  the 
early  work,  which  is  aligned  with  the  Tomahawk-to-task  allo¬ 
cation  carried  out  by  the  Tactical  Tomahawk  Weapon  Control 
System  (TTWCS)  Urioste  [34] . 

Weapon  assignment  is  part  of  a  large  class  of  opti¬ 
mization  models  known  as  resource  allocation  problems. 
Jaiswal  [19,  pp.  132-137,  and  references  therein]  divides  the 
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Figure  1.  Tactical  Tomahawk  Missile.  [Source:  Naval  Surface  Warfare  Center,  Dahlgren  Division,  2003]. 


complex  decision  process  of  planning  these  missions  into 
three  hierarchically  structured  classes:  (a)  weapon  mix  prob¬ 
lems,  which  determine  the  acquisition  of  a  weapons  mix  that 
inflicts  maximum  damage  to  a  set  of  targets,  subject  to  con¬ 
straints  including  cost,  manpower,  and  weapon  availability; 
(b)  weapon  deployment  problems,  which  assume  a  given 
weapon  mix  and  candidate  locations  where  these  weapons 
can  be  deployed,  and  which  seek  a  weapon  deployment  to 
maximize  damage  on  a  given  set  of  targets;  and  (c)  the 
weapon-to-target  allocation  problem,  which  is  “Given  the 
mix  of  weapons  deployed  at  specified  locations,  determine 
the  allocation  of  weapons  to  each  target  to  inflict  desired 
damage.”  The  resource  allocation  concepts  from  our  literature 
review  have  been  modified,  abstracted,  and  extended  in  order 
to  conform  to  the  task-to-firing  unit  allocation  addressed  in 
this  article. 

2.  THE  TOMAHAWK  MISSILE  AND  ITS 
TACTICAL  USE 

The  Tomahawk  entered  service  in  1983,  and  first  saw  com¬ 
bat  in  the  1 990-199 1  Gulf  War.  Tomahawks  are  now  available 
in  a  number  of  variants,  differing  in  the  type  of  munitions  car¬ 
ried,  guidance  system,  and  other  features.  Newer  versions 
have  improved  satellite  communications  that  allow  repro¬ 
gramming  the  missile  while  in  flight,  faster  launch  times,  and 
the  ability  to  loiter  in  an  area  of  emerging  targets.  Tomahawks 
use  both  inertial  and  GPS  navigation  systems,  as  well  as  ter¬ 
rain  contour  matching  radar  guidance  [13]  that  compares  a 
stored  map  reference  with  the  actual  terrain  to  determine  the 
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missile’s  position.  A  digital  optical  system  that  compares  a 
stored  image  of  the  target  with  the  observed  target  provides 
guidance  in  the  terminal  target  area. 

The  Tomahawk  has  proven  to  be  a  remarkably  flexible 
weapon.  In  1991,  hundreds  of  Tomahawks  were  launched 
against  land  targets  in  Iraq  during  Operation  Desert  Storm , 
principally  to  attack  the  Iraqi  air  defense  network.  In 
September  1995,  USS  Normandy  (CG  60)  fired  13 
Tomahawks  against  Bosnian  Serb  targets,  contributing 
to  the  subsequent  cease-fire  in  Bosnia  [3].  Tomahawks 
were  used  again  in  1996  during  Operation  Desert  Strike 
(Fig.  2). 


Figure  2.  Flying  her  battle  flag,  USS  Shiloh  (CG-67)  fires  a 
Tomahawk  land  attack  missile  in  Operation  Desert  Strike, 
3  September  1996. 
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initial  loadouts  on  two  AEGIS  VLS  half-modules 


Destroyer 

Tl 

Tl 

T2 

T3 

Cruiser 

Tl 

T2 

T2 

T3 

Loadoute  After  mmial  allocation  L  aado  uts  after  optimal  allocation 


Destroyer 

Tl 

Tl 

F 

Destroyer 

Tl 

72 

73 

Cruiser 

T2 

T2  \  Ti 

Cruiser 

Tl 

72 

T3 

Figure  3.  Impact  of  a  Tomahawk  missile  sequence  number  task-to-firing  unit  allocation  on  subsequent  strike  plan  capability. 


Equipping  submarines  with  the  highly  accurate  Tomahawk 
adds  a  new  conventional  land  attack  strike  warfare  capabil¬ 
ity.  The  employment  of  Tomahawks  enables  surface  ships  and 
submarines  to  hit  a  myriad  of  key  targets  such  as  command 
centers,  fortified  positions,  and  airfields. 

Before  leaving  port,  naval  ships  and  submarines  (referred 
to  as  the  “firing  units”  in  this  article)  are  loaded  with  Tom¬ 
ahawk  missiles.  The  fleet  command  determines  the  num¬ 
ber  of  Tomahawks  and  the  specific  variants  loaded  on  a 
firing  unit.  The  weapon  station’s  (port  facility)  and  firing 
unit’s  commander  decide  where  each  specific  Tomahawk  is 
loaded  in  the  launcher  or  torpedo  room.  For  AEGIS  cruisers 
and  destroyers  [20],  the  Naval  Surface  Warfare  Center  Port 
Hueneme  Division  designed  and  maintains  a  computer  pro¬ 
gram  that  suggests  the  exact  Tomahawk  loadout  positions. 
The  program  takes  into  account  the  capabilities  and  limi¬ 
tations  of  the  vertical  launch  system  (VLS)  on  the  AEGIS 
ships. 

A  Tomahawk  Missile  Sequence  Number  (MSN)  task-to- 
firing  unit  allocation  is  decided  at  the  fleet  command  level. 
This  allocates  each  Tomahawk  MSN  task  in  a  strike  plan 
to  a  particular  firing  unit;  a  missile-to-Tomahawk  MSN  task 
allocation  on  a  firing  unit  is  assigned  by  Tactical  Tomahawk 
Weapons  Control  System  on  each  firing  unit.  This  allocates 
a  specific  Tomahawk  loaded  on  the  firing  unit  to  each  par¬ 
ticular  Tomahawk  MSN  task  assigned  to  that  firing  unit;  this 
capability  has  been  documented  in  a  number  of  Naval  Surface 
Warfare  Center  Dahlgren  Division  (NSWCDD)  [27]  techni¬ 
cal  reports  (see,  e.g.,  Refs.  1, 14  and  15)  and  is  not  discussed 
further  in  this  article. 

For  a  given  set  of  Tomahawk  MSN  tasks  in  a  strike  plan, 
a  Tomahawk  Strike  Coordinator  (TSC),  usually  a  Lieutenant 
Commander  or  Commander,  develops  one,  or  more,  Tom¬ 
ahawk  MSN  task-to-firing  unit  allocation  alternatives  and 
presents  these  to  the  fleet  command  for  approval.  The  TSC 
must  consider  and  balance  many  factors  in  building  a  strike 
plan.  These  include  geographical  (e.g.,  distance  to  targets 
from  available  firing  units,  aim  points,  flight  route  coordina¬ 
tion),  logistical  (e.g.,  inventory  of  missiles  available,  techni¬ 
cal  conflicts  with  other  missiles,  preparation  time,  required 
time  on  target),  and  operational  (e.g.,  need  for  redundancy 
in  case  of  failure)  factors.  Simultaneously,  the  TSC  must 


consider  many  potentially  conflicting  objectives  and  prefer¬ 
ences,  such  as  [7]: 

•  Allocating  tasks  to  firing  units  so  that  each  firing  unit 
has  a  Tomahawk  variant  that  can  perform  the  task  allo¬ 
cated  in  a  cost-effective  manner.  For  example,  if  the 
target  can  be  struck  by  a  Tomahawk  variant  with  ter¬ 
rain  contour  matching  radar  guidance  capability,  we 
prefer  not  to  launch  an  “overqualified”  variant  which 
also  has  GPS  navigation  capability. 

•  Leaving  each  firing  unit  with  a  sufficiently  large  resid¬ 
ual  strike  capacity,  both  individually  and  as  part  of  a 
battle  group,  so  that  it  continues  to  be  effective  as 
a  naval  force  after  expending  the  assigned  missiles. 
Firing  unit  assignments  allocate  missiles  to  existing 
targets  while  platform  residual  capacity  preserves  the 
ability  to  address  potential  targets  (as  exemplified  in 
the  example  following  this  bulleted  list). 

•  Expending  missiles  first  from  firing  units  that  will 
soon  depart  the  theater  of  operations,  and  avoiding 
assignments  that  interfere  with  firing  units  engaged 
in  other  duties. 

The  TSC  has  a  very  difficult  job  that  must  often  be  per¬ 
formed  under  great  time  pressure,  and  in  conjunction  with 
other  duties.  As  an  illustration  of  the  complexity  of  the  prob¬ 
lem,  consider  a  simple  example  in  which  there  are  only  two 
firing  units  in  a  naval  strike  force:  an  AEGIS  destroyer  and 
an  AEGIS  cruiser.  On  both  AEGIS  platforms,  missiles  are 
carried  in  VLS  modules  that  can  contain  eight  cells,  each  of 
which  can  hold  a  Tomahawk.  Each  module  is  divided  into 
two  half-modules  of  four  cells.  At  most  one  Tomahawk  can 
be  powered  up  at  a  time  in  a  VLS  half-module.  For  this  exam¬ 
ple,  Fig.  3  shows  the  impact  of  a  Tomahawk  missile  sequence 
number  task- to- firing  unit  allocation  on  subsequent  strike 
plan  capability.  The  upper  part  of  the  figure  shows  hypo¬ 
thetical  initial  loadouts  in  a  VLS  half-module  on  the  AEGIS 
destroyer  and  an  AEGIS  cruiser.  The  AEGIS  destroyer  is 
loaded  with  four  Tomahawks  in  a  single  half-module  of  the 
firing  unit’s  VLS  (and  no  other  Tomahawks).  Specifically, 
there  are  two  variant  Tl  Tomahawks,  one  variant  T2  Tom¬ 
ahawk,  and  one  T3  in  the  half-module.  The  AEGIS  cruiser 
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is  loaded  with  four  Tomahawks  in  a  single  half-module  of 
the  firing  unit’s  VLS  (and  no  other  Tomahawks).  There  is 
one  variant  Tl  Tomahawk,  two  T2s,  and  one  T3  in  the 
half-module. 

The  TSC  creates  a  strike  plan  that  has  two  MSN  tasks. 
The  first  MSN  task  (MSN  1)  requires  the  launch  of  a  single 
missile,  and  can  be  performed  by  either  a  Tomahawk  variant 
Tl  or  a  Tomahawk  variant  T2;  but,  it  is  more  cost  effective 
to  use  a  T2  because  the  T 1  has  all  the  capabilities  of  a  T2  as 
well  as  additional  capabilities  that  are  not  needed  to  perform 
MSN  1.  Thus,  using  an  overqualified  Tl  would  waste  capa¬ 
bility  and  might  degrade  the  readiness  of  the  force  to  launch 
future  strikes. 

The  Tomahawk  selected  to  perform  MSN  1  must  be  pow¬ 
ered  up  and  launched  at  the  same  time  as  the  Tomahawk 
launch  time  for  the  other  MSN  task  (MSN  2).  This  second 
task  requires  the  launch  of  a  single  Tomahawk  variant  Tl. 
Because  both  MSN  tasks  have  the  same  missile  launch  time, 
they  conflict,  and  may  not  be  launched  from  the  same  plat¬ 
form.  Therefore,  one  task  must  be  allocated  to  the  AEGIS 
destroyer,  and  one  task  to  the  AEGIS  cruiser. 

Suppose  we  assign  MSN  1  to  the  AEGIS  destroyer  and 
MSN  2  to  the  AEGIS  cruiser.  Then,  onboard  the  AEGIS 
destroyer,  the  Tactical  Tomahawk  Weapon  Control  System 
(TTWCS)  missile-to-Tomahawk  MSN  task  allocation  would 
assign  the  Tomahawk  variant  T2  to  MSN  1  in  order  to  max¬ 
imize  the  future  launch  capability  of  the  AEGIS  destroyer. 
Onboard  the  AEGIS  cruiser,  the  Tomahawk  variant  T 1  would 
be  assigned  to  MSN  2  because  it  is  the  only  Tomahawk 
missile  that  can  perform  the  task.  Notice  that  as  a  result  of 
this  allocation  it  would  not  be  possible  to  perform  a  sub¬ 
sequent  strike  plan  requiring  two  Tls  to  be  launched  at 
the  same  time,  nor  a  strike  plan  that  might  be  performed 
by  launching  two,  more  cost-effective,  T2s  at  the  same 
time. 

However,  if  we  allocate  MSN  1  to  the  AEGIS  cruiser,  and 
MSN  2  to  the  AEGIS  destroyer,  then  onboard  the  AEGIS 
destroyer  the  TTWCS  missile-to-Tomahawk  MSN  task  allo¬ 
cation  would  assign  a  Tl  to  MSN  2;  and,  onboard  the  AEGIS 
cruiser  a  T2  would  be  assigned  to  MSN  1  because  the  Tom¬ 
ahawk  variant  T2  is  more  cost-effective  to  use  than  the 
Tomahawk  variant  Tl.  Notice  that  the  result  of  this  alloca¬ 
tion  allows  us  to  execute  a  subsequent  strike  plan  that  requires 
two  Tomahawk  variant  Tls  to  be  launched  at  the  same  time, 
or  a  subsequent  strike  plan  in  which  two  more  cost-effective 
missiles  can  be  launched  at  the  same  time. 

While  it  might  be  easy  to  find  the  preferred  Tomahawk 
MSN  task-to-firing  unit  allocation  in  this  simple  example  by 
inspection,  this  is  too  complex  to  do  when  multiple  Toma¬ 
hawk  MSN  tasks  with  conflicting  times  need  to  be  allocated 
to  multiple  firing  units  in  a  naval  strike  force.  As  indicated 
in  the  example  above,  the  missile-to-Tomahawk  MSN  task 
allocation  capabilities  that  were  (a)  initially  implemented  by 
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Naval  Surface  Warfare  Center,  Port  Hueneme  Division,  as 
a  stand-alone  capability,  (b)  subsequently  implemented  by 
NSWCDD  within  TTWCS,  and  (c)  later  formulated  indepen¬ 
dently  by  Kuykendall  [25],  all  reflect  that  the  missile-to-task 
allocation  on  board  a  firing  unit  has  to  be  performed  in  a 
manner  that  maximizes  future  launch  capability  of  the  firing 
unit. 

3.  PROBLEM  DEFINITION 

We  first  provide  some  definitions  and  rules  needed  to  define 
the  problem  and  formulate  a  model. 

3.1.  Conflicting  Tasks 

Two  Tomahawk  MSN  tasks  are  said  to  be  in  conflict  if  the 
preparation  and  launch  time  epochs  for  each  missile  allocated 
to  perform  the  Tomahawk  MSN  tasks  overlap.  Time  conflict 
task  sets  are  constructed  based  on  the  preparation  and  launch 
time  epochs  of  the  Tomahawk  MSN  tasks,  and  the  various  fir¬ 
ing  unit  capabilities  and  restrictions.  A  Tomahawk  MSN  task 
may  be  in  more  than  one  time  conflict  task  set.  Constraints 
are  generated  to  reflect  each  firing  unit’s  capability  to  handle 
tasks  with  conflicting  times. 

3.2.  Firing  Unit  Characteristics 

There  are  several  types  of  surface  ship  and  submarine  firing 
units.  Each  type  of  firing  unit  has  particular  weapons  storage 
and  launch  characteristics  that  impose  different  constraints. 
The  surface  ships  store  their  Tomahawks  in  VLSs.  Each  ship 
may  hold  a  different  number  of  missiles,  and  the  cell-location 
of  the  missiles  may  differ  from  ship  to  ship.  Depending  on  the 
class  of  submarine,  Tomahawks  can  be  stored  in  the  torpedo 
room  or  in  a  vertical  launch  canister,  with  capacity  varying 
by  submarine  class. 

Allocating  a  Tomahawk  MSN  task  to  a  firing  unit  that 
is  already  engaged  in  some  other  activity  incurs  a  so-called 
employment  penalty;  thus,  an  otherwise  engaged  firing  unit 
is  used  to  perform  such  a  task  only  if  necessary. 

Each  firing  unit  is  characterized  as  either  an  expend  or  a 
non-expend  unit.  A  firing  unit  might  be  designated  as  expend 
if  it  is  scheduled  to  leave  the  theater  of  operations  soon,  and 
thus  it  is  advantageous  to  the  remaining  force  in  theater  if  the 
departing  firing  unit  launches  its  weapons.  For  non-expend 
firing  units,  we  allocate  Tomahawk  MSN  tasks  to  maximize 
the  ability  to  perform  follow-on  strikes. 

3.3.  Tasks  and  Task-Parts 

A  Tomahawk  MSN  “task”  is  made  up  of  one  or  more 
missiles  that  are  to  be  prepared  for  launch  to  fly  a  partic¬ 
ular  over-land  route  and  attack  a  specified  aimpoint.  A  task 
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can  have  one  or  more  task-parts  (defined  below),  each  of 
which  requires  at  least  one  missile.  “Strike  packages”  are 
sets  of  tasks  or  missions  created  by  the  TSC.  For  exam¬ 
ple,  a  strike  package  may  contain  the  following  information: 
“Task  1,  requires  one  primary  part,  one  ready-spare  part,  and 
one  backup  part,  Tomahawk  variant  Tl ,  with  earliest/latest 
time  to  launch  00:35/00:40;  Task  2,  requires  one  primary 
part  and  one  backup  part,  Tomahawk  variant  T2,  with  ear¬ 
liest/latest  time  to  launch  00:37/00:42;  etc.”  After  the  strike 
packages  are  created,  the  TSC  sends  them  (via  the  designated 
communication  system)  to  the  firing  units.  The  firing  units 
return  exception  reports  detailing  missions  they  cannot  com¬ 
plete  as  specified.  Note  that  in  our  system,  the  allocation  of 
successive  strike  packages  is  dynamically  incorporated  into 
the  optimization.  That  is,  the  execution  sequence  freezes  the 
assets  selected  for  previous  allocations  and  proceeds  with 
subsequent  packages  from  the  remaining  assets. 

There  are  five  types  of  task-parts: 

•  A  primary  task-part  requires  one  or  more  primary 
missiles  that  are  to  be  launched  to  fly  a  particular 
over-land  route  and  are  to  attack  a  specified  aimpoint. 

•  A  ready-spare  task-part  is  associated  with  a  specific 
primary  task-part,  and  requires  one  or  more  ready- 
spare  missiles  (located  on  the  same  firing  units  as 
the  task’s  primary  missiles).  The  ready-spare  mis¬ 
siles  provide  redundancy  for  the  primary  missiles, 
in  case  some  of  the  primary  missiles  do  not  launch 
successfully. 

•  A  backup  task-part  is  associated  with  a  specific  pri¬ 
mary  task-part,  and  requires  one  or  more  backup 
missiles  located  on  firing  units  other  than  those  that 
have  the  task’s  primary  missiles.  The  backup  missiles 
provide  a  second-level  redundancy  for  the  primary 
missiles,  in  case  the  required  total  number  of  primary 
and  ready-spare  missiles  is  not  launched  successfully 
from  primary  firing  units. 

For  completeness,  we  describe  two  other  types  of  task- 
parts  that  are  set  aside  at  the  outset,  that  is,  they  are  not 
available  in  the  allocation  problem  modeled  in  this  article. 

•  A  call-for-fire  task-part  requires  one  or  more  missiles 
be  reserved  for  launching  in  support  of  friendly  force 
shore  operations  during  a  specified  time  interval. 

•  A  pooled  task-part  designates  one  or  more  missiles 
to  be  ready  to  be  used  in  support  of  last-minute  tasks 
during  a  specified  time  interval.  If  a  last-minute  need 
arises  for  a  firing  unit  with  pooled  missiles  to  perform 
a  specific  primary  task,  then  the  appropriate  number 
of  pooled  missiles  can  be  selected  for  launch  from  the 
firing  unit  to  attack  the  specified  aimpoint. 


A  number  of  constraints  guide  the  assignment  of  task- 
parts: 

•  Tasks  are  classified  as  separable  or  nonseparable :  if 
a  task  is  separable  (nonseparable)  then  all  the  missile 
requirements  of  a  task-part  associated  with  the  task 
can  (cannot)  be  allocated  to  different  firing  units. 

•  Some  firing  units  are  designated  as  spread ,  and,  if  pos¬ 
sible,  the  primary  and  back-up  task-part  assignments 
should  be  allocated  among  as  many  of  these  firing 
units  as  possible.  Note:  A  related  concept  is  known 
as  “task-part  spreading.”  For  example,  a  primary  task- 
part  that  is  designated  as  spread  should  be  allocated 
among  different  firing  units.  This  article  incorporates 
spread  firing  units,  but  not  spread  task-parts. 

•  A  ready-spare  task-part  must  be  assigned  to  a  firing 
unit  that  is  assigned  an  associated  primary  task-part. 
Also,  in  the  allocation  of  ready- spare  and  primary 
task-parts'  for  a  specific  separable  task,  for  each  fir¬ 
ing  unit,  the  required  number  of  ready-spare  missiles 
needed  for  the  ready-spare  task-part  should  be  no 
more  than  the  number  of  primary  missiles  needed  for 
the  primary  task-part. 

•  A  backup  task-part  must  not  be  assigned  to  a  firing 
unit  that  is  assigned  an  associated  primary  task-part. 
Some  firing  units,  especially  submarines,  may  be 
excluded  from  being  assigned  backup  task-parts'  by 
the  TSC. 

•  A  task  is  considered  “allocated”  only  when  all  the 
missiles  required  for  each  of  its  task-parts  have  been 
assigned.  Partial  allocation  of  a  task  is  not  allowed 
in  practice  and  our  formulation  explicitly  avoids  that 
possibility. 

3.4.  Launch  Areas 

Each  firing  unit  is  located  in  a  specific  launch  area,  and 
can  only  be  assigned  to  strike  targets  that  are  geographically 
feasible  from  that  area.  Relocation  of  the  firing  unit  to  a  dif¬ 
ferent  launch  area  is  beyond  our  planning  scope,  and  is  not 
considered. 

3.5.  Missile-to-Mission  Matching 

There  are  multiple  Tomahawk  missile  variants,  and  each 
one  is  designated  by  its  Expanded  Missile  Identification 
(XMID).  For  the  missile-to-task  allocation  that  takes  place 
on  the  individual  firing  units,  the  Tomahawks  that  are  capa¬ 
ble  of  performing  each  task  are  identified,  prioritized,  and 
allocated  based  on  a  number  of  factors.  Although  the  TSC 
does  not  have  information  concerning  all  of  the  factors,  he 
has  knowledge  concerning  the  Tomahawk  missile  loadout  on 
each  firing  unit,  and  the  cost  effectiveness  of  using  a  particu¬ 
lar  Tomahawk  variant  (i.e.,  XMID)  to  perform  a  task.  For  any 
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task,  a  Missile  Mission  Matching  (or  M3)  list  of  XMIDs  that 
are  able  to  perform  the  task  is  constructed  and  prioritized  from 
the  most  to  the  least  cost-effective  Tomahawk  variant  to  per¬ 
form  the  task.  The  task-to-firing  unit  allocation  should  assign 
task-parts  to  firing  units  in  order  to  maximize  the  potential 
follow-on  launch  capability  of  the  force  by  allowing  the  fir¬ 
ing  units  to  perform  Tomahawk-to-task  allocations  that  use 
the  most  cost  effective  Tomahawks. 

3.6.  Tomahawk  Strike  Coordinator  Overrides 

While  our  intent  here  is  to  develop  an  automated  deci¬ 
sion  aid  to  allocate  task-parts  to  firing  units,  there  are  cases 
in  which  the  suggested  allocation  should  not  be  applied  and 
other  cases  in  which  the  TSC  may  wish  to  override  the  task 
part  to  firing  unit  allocation,  e.g.: 

•  Tasks  can  be  divided  into  two  groups:  execute  tasks 
and  launch-plan  tasks.  “Execute”  tasks  are  tasks 
whose  parts  have  already  been  allocated  to  firing 
units,  and  the  firing  units’  commanders  have  con¬ 
firmed  that  Tomahawks  can  be  allocated  to  the  task 
parts.  Within  the  automated  decision  aid  to  allocate 
task-parts  to  firing  units,  any  prior  execute  task  allo¬ 
cations  to  firing  units  should  be  accounted  for,  but 
not  changed.  “Launch-plan”  tasks  are  those  whose 
task  parts  have  not  yet  been  allocated  to  firing  units. 
The  job  of  the  automated  decision  aid  is  to  allocate 
launch-plan  task-parts  to  firing  units,  while  taking 
prior  execute  tasking  into  account. 

•  The  TSC  can  allocate  a  particular  task-part,  or  a  por¬ 
tion  of  a  task-part,  to  a  specific  firing  unit,  or  even 
designate  launch  of  a  specific  missile.  Such  general, 
ad  hoc  constraints  offering  the  TSC  complete  control 
over  predesignation  are  a  key  to  gaining  and  keeping 
the  confidence  and  trust  of  the  TSC.  In  our  experi¬ 
ence,  planners  distrust  “black  box”  tools  over  which 
they  have  no  control,  and  find  ways  to  avoid  using 
such  tools. 

4.  MATHEMATICAL  MODELING 
4.1.  Objective  Functions 

We  have  developed  objective  functions  by  interviewing 
TSCs,  their  commanding  officers  and  senior  naval  officials. 
The  objective  functions  listed  below  are  for  a  streamlined 
model  formulation  presented  in  the  sections  that  follow  that 
deals  only  with  the  allocation  of  Tomahawk  tasks  to  AEGIS 
ships.  The  notation  is  defined  in  detail  in  the  Appendix,  but 
the  definitions  of  various  terms  are  repeated  here. 

OBJECTIVE  1  (Number  of  unassigned  task-parts;  to  be 
minimized):  This  objective  expresses  the  number  of  task- 
parts  not  allocated  to  firing  units.  While  partial  allocation 
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of  a  task  is  not  allowed,  this  objective  could  be  modified  to 
reflect  the  importance  and  relative  priority  of  the  various  tasks 
and  task-parts. 

Objective  k  =  1  :  min  z\  =  Utpm, 

teT,peP, 

meMtp 

where  Utpm  equals  1  if  the  m-th  missile  requirement  for 
task-part  p  in  task  t  is  unmet,  and  0  otherwise. 

OBJECTIVE  2  (Use  of  firing  units  currently  engaged  in 
other  tasking;  to  be  minimized):  This  objective  penalizes 
distraction  of  firing  units  engaged  in  other  missions. 

Objective  k  —  2  :  min  zi  =  psYs, 

seS 

where  Vs  equals  1  if  ship  s  is  used,  and  0  otherwise  and  where 
ps  is  the  penalty  for  using  ship  s. 

OBJECTIVE  3  (Number  of  primary  task-parts  allocated 
to  designated  “expend”  firing  units;  to  be  maximized):  This 
objective  rewards  firing  units  regarded  as  “expend”  if  they 
are  assigned  primary  task-parts  (not  redundancies). 

Objective  k  —  3  :  max  z3  =  Ystpm 

t  dr  ,s  G>Sf  n>s^ , 
peP,meMtp 

where  Ystpm  equals  1  if  ship  s  is  assigned  to  perform  the  m-th 
missile  requirement  for  task-part  p  in  task  t,  and  0  otherwise. 

OBJECTIVE  4  (Deviations  in  residual  missile  inventory 
across  firing  units;  to  be  minimized):  This  objective  seeks  to 
level  residual  missile  inventories  across  firing  units  (except 
for  expend  ships)  after  primary  task-parts  are  assigned.  The 
rationale  for  this  objective  is  to  avoid  weakening  non-expend 
ships  for  potential  future  engagements. 

Objective/:  =  4  :  min  za  =  ^  RDV, 

seS\SE 

where  RDS  is  the  difference  between  the  residual  number  of 
missiles  on  ship  s  and  the  average  residual  number  of  missiles 
across  all  “non-expend”  ships. 

OBJECTIVE  5  (Number  of  “spread”  firing  units  used 
for  primary  task-part  assignments;  to  be  maximized):  This 
objective  counts  the  number  of  firing  units  included  in  our 
predesignation.  Maximizing  this  number  spreads  the  alloca¬ 
tion  of  primary  task-parts,  minimizes  the  impact  of  losing 
one  or  several  ships,  and  encourages  participation  by  many 
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firing  units.  Some  ships  can  be  marked  to  be  “preferred”  for 
primary  task-parts  by  the  TSC. 

Objective  k  =  5  :  max  Z5  =  ^  PR^, 

seSPR 

where  PRS  equals  1  if  any  primary  task-part  is  allocated  to 
ship  s  (only  for  s  e  SFR),  and  0  otherwise.  SFR  c  S  is  a  sub¬ 
set  of  AEGIS  ship  firing  units  to  which  we  wish  to  allocate 
primary  task-parts. 

OBJECTIVE  6  (Number  of  “spread”  firing  units  used 
for  backup  task-part  assignments;  to  be  maximized):  This 
objective  is  the  same  as  Objective  5,  but  applies  to  backup 
task-parts. 

Objective  k  =  6  :  max  Z6  =  BU^, 

seSBU 

where  BU*  equals  1  if  any  back-up  task-part  is  allocated  to 
ship  s  (only  for  s  e  SBU),  and  0  otherwise.  SBlJ  c.  S  is  a  sub¬ 
set  of  AEGIS  ship  firing  units  to  which  we  wish  to  allocate 
backup  task-parts. 

OBJECTIVE  7  (Use  of  non-cost-effective  missiles;  to  be 
minimized):  This  objective  assesses  the  cost  effectiveness 
of  the  missiles  that  may  be  allocated  to  tasks  on  the  firing 
units,  and  minimizing  the  use  of  non-cost-effective  missiles 
preserves  potential  follow-on  launch  capability. 

Objective  k  —  1:  min  zi  =  ^  m3wpXshwtptn, 

t  G  T  ,s  ,h€zHs , 
weWsht,peP,meMtp 

where  ,  is  the  penalty  for  allocating  weapon  w  to  part  p 
in  task  t  (based  on  M3  position  of  weapon  w  for  task  t ,  and  the 
weight  of  part  p  to  match  the  M3  list).  Also,  Xshwtpm  equals  1 
if  ship  s  is  assigned  to  perform  the  m-th  missile  requirement 
for  task-part  p  in  task  t  (tentatively  using  a  weapon  of  type  w 
in  half-module  h ),  and  0  otherwise. 

OBJECTIVE  8  (Relative  value  of  residual  salvo;  to  be 
maximized):  This  objective  assesses  the  potential  missile 
loadout  on  the  firing  units  after  Tomahawks  are  launched  to 
perform  primary  task  parts,  in  order  to  maximize  the  potential 
follow-on  salvo,  that  is,  launch  capability.  (The  formulation 
presented  below  only  considers  the  “most  valuable”  missile 
type  in  a  launcher  half  module.  The  objective  can  be  modified 
to  reflect  all  the  missile  types  in  a  launcher  half  module.) 

Objective  k  =  8  :  max  z8  =  ^  vwDshw , 

S  £  S,h  £  Hs ,  W  £  Wsh 

where  vw  is  the  relative  value  of  weapon  w  for  residual  salvo 
computation. 


Also,  if  based  on  the  information  held  by  the  TSC,  one 
or  more  weapons  of  type  w  would  remain  on  ship  s  in  half¬ 
module  h  after  firing  all  primary  task-part  missiles,  at  most  1 
Dshw  for  that  (s,h)  is  allowed  to  be  equal  to  1,  and  must  be  0 
otherwise. 

4.2.  Model  Formulation 

We  present  a  simplified  model  formulation  that  covers 
most  problem  features.  In  this  streamlined  formulation,  we 
have  suppressed  geographic  proximity,  substituting  instead 
lists  of  geographically  feasible  assignments.  Here,  we  con¬ 
sider  only  AEGIS  ships  as  firing  units.  For  details  on  mod¬ 
eling  with  SSN  submarines,  see  Ref.  9.  For  details  on  SSGN 
modeling,  see  Ref.  8.  The  inclusion  of  submarines  does  not 
affect  the  solvability  of  the  problem  any  more  than  adding 
new  surface  ships.  Zumwalt-class  ships  [33]  have  not  yet  been 
inserted  in  the  formulation.  (For  specifications  of  various 
combatant  ship  classes,  see  Ref.  21.)  Only  primary,  ready- 
spare,  and  backup  task-parts  are  discussed  in  the  following 
formulation.  A  structured  list  of  variable  and  parameter  defi¬ 
nitions  and  a  formal  definition  of  constraints  are  provided  in 
the  Appendix. 

4.3.  Solving  the  Multiobjective  Optimization  Problem 

The  model  that  we  have  described  is  a  mixed-integer  lin¬ 
ear  program  with  multiple  objective  functions.  In  general,  a 
multiobjective  optimization  problem  may  be  written  as  (e.g., 
Ref.  11,  p.  16): 

“minimize” F[zi(x),  Z2M, . . .  ,Z|*|(*)] 
subject  to  x  e  0 

where  x  =  {x\,  *2, . . . ,  xn}  is  a  vector  of  decision  variables, 
and  zu-.-,Z\k\  are  objective  functions  of  these  decision 
variables,  0  is  the  set  of  feasible  decisions  allowed  by  the 
constraints,  and  F  is  some  vector  or  scalar  function  of  the 
objective  function  values.  By  necessity,  the  term  “minimize” 
may  have  several  meanings  and  the  choice  among  them  is  sub¬ 
jective.  If  there  is  a  value  of  x  that  simultaneously  optimizes 
each  ZkM,  then  the  multiobjective  problem  would  be  clear 
and  simple,  but  such  an  ideal  solution  cannot  be  guaranteed, 
and  is  rare  in  practice. 

For  convenience,  we  assume  that  any  objective  that  is  to 
be  maximized  in  the  original  formulation  is  redefined  as  its 
negation,  so  that  all  objectives  are  to  be  minimized. 

4.3.1.  The  Value  Function 

If  F  is  a  well-defined,  scalar- valued  function  that  all  inter¬ 
ested  parties  agree  is  an  accurate  representation  of  their 
preferences  and  priorities,  then  the  multiobjective  problem 
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reduces  to  single-objective  optimization,  for  which  standard 
methods  can  be  employed.  Such  F  functions  are  called  value 
functions  in  a  deterministic  setting  (which  is  the  case  here)  or 
they  are  called  utility  functions  in  a  stochastic  setting.  (See 
Keeney  and  Raiffa  [22,  pp.  82-116  and  131-218,  respec¬ 
tively]  for  development  of  these  functions  and  Harrison  and 
Rosenthal  [16]  for  their  use  in  optimization  models.) 

4.3.2.  Linear  Value  Function 

The  most  commonly  used  value  function  is  a  weighted 
sum  of  the  objective  functions  that  has  the  advantage  of 
being  intuitively  appealing  to  many  people.  This  also  has 
well-documented  drawbacks.  The  problem  becomes 
intractable  if  a  nonlinear  weighting  is  used.  Despite  its  poten¬ 
tial  weakness,  a  linear  value  function  (that  is,  weighted  sums) 
is  widely  used  and  is  often  an  advisable  approach,  especially 
when  under  time  pressure.  There  is  no  way  to  avoid  some 
subjective  judgment  on  empirical  observations  of  prelimi¬ 
nary  results  when  setting  the  weights.  We  look  at  the  need 
for  subjective  judgment  in  multiobjective  optimization  (and 
other  areas  of  OR)  as  an  important  distinguishing  feature 
of  OR  versus,  say,  physics,  because  “we  are  engaged  in  the 
modeling  of  systems  that  involve  people,  not  particles.”1 

The  formulation  of  the  linear- value-function  or  weighted- 
sum  is: 

(W)  :  Z*  =  min  za  =  '22  MkZk(x) 

keK 

S.t.  X  €  6 

where  cok  >  0  is  the  weight  given  to  objective  k. 

4.3.3.  Pareto  Optimality 

For  a  multiobjective  problem,  it  is  always  preferred  to 
choose  a  Pareto-optimal  solution  (e.g.,  Ref.  [11],  p.  24),  that 
is,  one  from  which  it  is  impossible  to  find  another  point  that 
improves  the  attainment  of  some  objective(s)  without  simul¬ 
taneously  diminishing  the  attainment  of  some  other  objective. 
There  are  usually  many  Pareto-optimal  points,  perhaps  an 
infinite  number.  Where  \K\  =  2,  the  Pareto-optimal  set  is 
a  curve,  the  efficient  frontier,  which  is  easy  to  generate  and 
extremely  useful  for  presenting  alternatives  to  decision  mak¬ 
ers.  The  case  of  |  K  \  >  2  is  more  difficult  and  has  been  treated 
by  Zeleny  [37],  Rosenthal  [28],  Steuer  [31],  and  White  [35]. 
One  fortunate  result  in  the  theory  of  multiobjective  opti¬ 
mization  is  Soland’s  Theorem  [30],  from  which  it  follows 
that  solving  problem  (W)  yields  a  Pareto-optimal  solution, 
although  not  all  efficient  points  may  be  derived  as  the  solution 
to  a  problem  like  (W). 

1  This  observation  is  classic  Rick  Rosenthal. 
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4.3.4.  Hierarchical  Optimization 

Another  intuitively  simple  idea  is  prioritizing  the  multiple 
objective  functions  and  then  optimizing  them,  one  at  a  time, 
in  priority  order.  Attainment  of  optimal  (or  near-optimal)  val¬ 
ues  of  higher  priority  objectives  becomes  a  constraint  when 
optimizing  the  lower  priority  ones.  This  method  requires 
solving  the  following  sequence  of  increasingly  restricted 
single-objective  problems  (H^),  for  k  =  1, . . . ,  \K\: 

(H*):  4=  min  Zk(x) 

S.t.  X  €  6 

zv(x)  =  4,  V*'  <  k 

One  weakness  of  this  hierarchical  optimization  (also  known 
as  pre-emptive  priorities)  is  that  the  optimization  at  each  step 
depends  upon  the  existence  of  multiple  optima  in  the  previous 
step.  Without  multiple  optima  in  step  k ,  the  method  con¬ 
verges  before  considering  objectives  &  +l,...,|K|.  Another 
valid  criticism  of  hierarchical  optimization  is  that  it  may 
bar  consideration  of  some  sensible  opportunities  to  trade 
off  a  small  decrement  in  a  high-priority  objective  in  return 
for  a  large  improvement  in  a  low-priority  one.  These  issues 
can  be  addressed  by  relaxing  the  higher-priority  objective’s 
aspirations,  replacing  each  equality  constraint  in  by  an 
inequality  bounded  by  some  given  fraction  of  the  optimal 
value  z*k. 

If  the  weights  are  chosen  such  that  co\  >>  cl>2  >>  •  •  •  >  > 
co\k\  ,  then  the  weighted- sum  objective  yields  the  same  solu¬ 
tion  as  hierarchical  optimization.  In  practice,  a  large  disparity 
in  the  magnitude  of  the  weights  and  objective  coefficients 
can  be  computationally  problematic.  This  is  because  in  order 
for  the  low-priority  objectives’  weights  to  have  any  influ¬ 
ence  on  the  solution,  the  optimality  tolerance  needs  to  be 
set  very  close  to  zero,  which  in  turn  may  cause  the  lower- 
priority  objectives’  contributions  to  the  value  function  to  be 
ignored  by  the  optimization  numerics.  Normalization  and 
scaling  across  all  objectives  should  be  implemented  to  help 
reduce  this  risk. 

The  (W)  and  (H)  formulations  above  can  be  solved  as  pre¬ 
sented,  or  approximated  by  a  problem  cascade  approach,  as 
in  Brown  et  al.  [6,  p.  341].  This  strategy  is  also  known  as 
fix-and-relax  [10, 12].  This  technique  relaxes  the  integrality 
of  most  of  the  decision  variables  to  make  the  problem  eas¬ 
ier.  For  example,  only  those  variables  related  to  a  subset  of 
tasks  are  forced  to  assume  integer  values.  Then,  these  vari¬ 
able  values  are  fixed  and  a  new  subset  of  variable  values  is 
forced  to  be  integral,  until  all  of  them  are  fixed.  Algorithmic 
details  are  provided  by  Demir  [9].  We  find  this  approach  may 
help  reduce  computational  time  for  cases  in  which  (W)  or  (H) 
take  a  long  time  to  converge,  without  much  degradation  in 
the  achieved  objectives.  However,  it  also  increases  the  com¬ 
putational  time  in  easier  instances  due  to  the  multiple  models 
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that  need  to  be  solved.  As  our  computational  experiments  are 
restricted  by  a  short  time  limit,  we  do  not  discuss  the  cascade 
approach. 

4.3.5.  Heuristic  Optimization 

In  addition  to  the  conventional  multiobjective  optimization 
methods  of  solving  (W)  and  (H^),  we  have  also  developed 
a  specialized,  custom  heuristic  for  finding  good  solutions 
quickly.  The  heuristic  uses  sequential  insertion:  Queues  of 
task-parts,  firing  units,  and  missiles  are  maintained,  with 
priorities  in  the  queues  determined  by  a  battery  of  problem- 
specific  criteria.  As  assignments  are  made  from  the  queues, 
the  queues  are  replenished  and  reordered.  Occasionally,  we 
invoke  a  restart  to  give  the  heuristic  an  opportunity  to  improve 
the  solution  with  large  changes,  rather  than  with  small  incre¬ 
ments.  Details  of  this  heuristic  appear  in  Salmeron  [29]  and 
Arnold  [2],  but  initial  steps  in  exploring  heuristics  for  this 
problem  are  due  to  Kirk  [23]  and  Hodge  [17]. 

Over  the  period  in  which  they  have  been  involved  in  mod¬ 
eling  Tomahawk  missile  assignment,  the  Naval  Postgraduate 
School  team  of  students  and  faculty  has  worked  closely  with 
experienced  operations  officers.  This  close  contact  with  oper¬ 
ational  practice  has  allowed  the  team  to  gain  domain- specific 
knowledge  that  allows  the  resolution  of  various  modeling 
issues.  For  example,  Kubu  [24]  systematically  studies  situ¬ 
ations  in  which  the  heuristic  is  unable  to  assign  all  tasks  to 
firing  units.  He  introduces  a  diagnostic  module  that  provides 
the  TSC  with  likely  reasons  as  to  why  a  task  cannot  be  cov¬ 
ered.  He  also  provides  a  prescriptive  module  for  suggesting 
changes  to  the  inputs  so  that  all  tasks  can  be  allocated.  Kubu 
collects  domain-specific  knowledge  about  our  problem  and 
assembles  a  catalog  of  suggested  adjustments  to  diagnose  and 
eliminate  infeasibilities.  The  most  effective  tool  is  an  adjust¬ 
ment  of  time  windows  performed  by  sliding  the  earliest  or 
latest  launch  times  in  order  to  avoid  task  conflicts. 

5.  TEST  PROBLEMS  AND  COMPUTATIONAL 
EXPERIENCE 

5.1.  Test  Problems 

We  have  developed  a  dataset  of  35  test  problems.  We  refer 
to  the  first  22  test  problems  as  lab  problems,  and  the  other  13 
as  exercise  problems.  The  lab  problems  are  complete  inven¬ 
tions  [2],  based  upon  war-fighter  and  engineering  judgment, 
but  without  empirical  basis.  In  contrast,  the  exercise  prob¬ 
lems  are  based  entirely  on  carefully  recorded  observations  of 
naval  exercises  at  sea  attended  by  Wingeart  [36].  Wingeart 
collected  both  inputs  and  outputs  of  all  the  tasks  to  firing  unit 
allocations  completed  in  two  major  U.S.  Navy  Third  Fleet 
exercises  known  as  Rimpac  and  Texas -Thunder.  His  analysis 


Table  1.  Predesignation  optimization  problem  characteristics  for 
lab  and  fleet  exercise  test  cases. 


No.  of 

No.  of 

No.  of 

Test  cases 

firing  units 

tasks 

tasks-parts 

22  Lab  problems  [2] 

4-8 

13-26 

104-232 

13  Exercise  problems  [36] 

4-8 

18-100 

28-128 

includes  comparisons  of  what  would  have  happened  had  our 
model  been  used,  versus  the  actual  exercise  results. 

Our  runs  were  made  on  a  2  GHz  laptop  with  2  GB  of 
RAM.  Exact  models  (H)  and  (W)  were  implemented  in 
GAMS  [5]  and  solved  with  CPLEX  9.0  [18].  The  heuristic 
was  implemented  in  Java  [32]. 

5.2.  Results  for  Lab  Problems 

Table  1  shows  problem  characteristics  for  both  lab  [2]  and 
Third  Fleet  exercises  [36]  test  cases.  Note  that  “x-y”  indicates 
a  range  from  a  minimum  of  x  to  a  maximum  of  y.  Tables  2 
and  3  show  detailed  results  for  the  22  lab  problems  and  the 
13  exercise  problems,  respectively. 

By  column,  Table  2  shows  the  case  name,  the  optimal  solu¬ 
tion  to  each  of  the  eight  goals,  and  the  computational  time  for 
each  of  the  three  methods:  Heuristic,  (W)  and  (H)  models, 
respectively.  For  all  of  these  test  cases,  the  exact  optimization 
solutions  (W)  and  (H)  are  identical.  In  some  cases,  the  heuris¬ 
tic  solution,  which  is  displayed  by  an  underlined  alternative 
value,  is  of  inferior  quality.  For  example,  all  three  methods 
achieve  the  same  solutions  for  cases  “1”  and  “1A,”  but  for 
case  “IB”  the  heuristic  solution  is  3  for  goal  k  —  6,  when 
the  optimal  objective  should  be  4.  Of  course,  the  heuristic 
solution  might  be  superior  for  goals  k  +  1, . . . ,  8  (see,  e.g., 
case  “2”)  but  that  can  only  happen  if  it  is  of  inferior  quality 
for  a  hierarchically  superior  goal  k.  Overall,  heuristic  solu¬ 
tions  are  of  acceptable  quality.  In  the  22  lab  problems,  the 
heuristic  never  fails  to  optimize  the  first  goal,  only  fails  to 
optimize  the  second  goal  in  one  of  the  problems  (case  “3”), 
and  the  third  goal  in  another  problem  (case  “6”).  Goal  k  =  4, 
“leveling  missiles  left  on  firing  units  by  minimizing  absolute 
deviations”  is  apparently  the  most  difficult  goal,  since  the 
heuristic  misses  the  optimal  solution  in  10  of  the  remain¬ 
ing  lab  cases.  Indeed,  most  of  the  time  consumed  by  the  (H) 
method  in  solving  the  eight  problems  (H^),  k  —  1, . . . ,  8,  is 
spent  in  solving  (H4). 

5.3.  Results  for  Exercise  Problems 

In  Table  3,  “R-Sft”  denotes  the  ft-th  strike  scenario  for 
the  Rimpac  exercise,  and  “T-S^-m”  denotes  the  m-th  salvo 
within  the  n-ih  strike  scenario  for  the  Texas  exercise.  (Rim¬ 
pac  strikes  have  one  unique  salvo  each.)  Strikes  are  planned 
independently  in  that  they  may  involve  different  groups  of 
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Table  2.  Result  comparison  for  lab  problems:  Heuristic  objective  value  (underlined,  if  different)  and  (W)-(H)  identical  objective  values. 


Test 

case 

Goals 

Solution  time 
(s):  Heuristic, 
(W),(H) 

#  1  (Min): 
Unassigned 
task  parts 

#  2  (Min): 

Penalty 
firing  units 

#  3  (Max): 
Expend 
units 

#  4  (Min): 
Deviation 
level 

#  5  (Max): 
Spread 
primary 

#  6  (Max): 
Spread 
back-up 

#  7  (Min): 
M3  list 
position 

#  8  (Max): 
Residual 
salvo 

1 

0 

0 

0 

34 

3 

0 

92 

160 

2,  1,  3 

1A 

0 

0 

10 

8.7 

0 

0 

92 

160 

2,  1,3 

IB 

0 

0 

10 

8.7 

0 

3,4 

92 

160 

2,  1,3 

2 

0 

0 

0 

27,  23 

2,0 

0 

221,  184 

162,  157 

3,  69,  4 

2A 

0 

1 

0 

44.5,  42.5 

4 

0 

196,  186 

141,  146 

4,  1,3 

3 

0 

1,0 

0 

1.5,  15.5 

3 

0 

368 

130,  140 

5,  8,  14 

4 

0 

1 

0 

16,0 

4 

0 

508,  483 

107,  127 

10,  12,  52 

5 

0 

0 

0 

25.4,  11.4 

6,4 

0 

599,  558 

252,  268 

10,  7,  32 

5A 

0 

0 

73 

40 

4 

0 

596 

227,  239 

10,  10,  42 

5B 

0 

1 

0 

35.1,  20.6 

5 

0 

629,  586 

240,  270 

16,  54,  329 

6 

0 

0 

69,78 

24,  13.5 

6,5 

0 

902,  882 

201,  240 

7,  42,  172 

7 

0 

0 

0 

36 

2,3 

0 

92 

160,  152 

2,  1,3 

7A 

0 

0 

5 

34.7,  28 

0 

0 

104,  92 

152 

5,  1,2 

7B 

0 

0 

10 

8.7 

0 

3 

92 

150 

2,  1,2 

8 

0 

0 

0 

36 

0 

0 

80 

162 

2,  1,2 

8A 

0 

0 

0 

36 

0 

0 

81,  80 

162,  153 

2,  1,2 

9 

0 

0 

0 

34 

3 

0 

92 

160 

3,  1,4 

9A 

0 

0 

0 

40,  36 

0 

3,2 

92 

160,  150 

3,  1,  1 

10 

0 

0 

0 

39,  35 

3 

0 

157, 149 

147,  151 

3,  1,4 

10A 

0 

0 

0 

37,  45 

3 

0 

167,  153 

147,  140 

3,  1,3 

11 

0 

0 

0 

34 

1,3 

0 

58,92 

162,  160 

1,  1,5 

12 

0 

0 

0 

43,  16.5 

7 

0 

654,  665 

256,  271 

12,  46,  206 

firing  units,  or  firing  units  may  reload  after  the  strike.  On 
the  other  hand,  salvos  within  a  strike  (usually  arranged  by 
earliest  launch  time)  are  planned  so  that  a  follow-on  salvo 
can  only  make  use  of  the  remaining  Tomahawk  inventory 


from  previous  salvos.  In  this  exercise,  the  fleet  assumes  a 
“best-case  strike  planning:”  only  primary  task-parts  are  actu¬ 
ally  fired,  whereas  redundant  allocations  become  available 
for  subsequent  salvos. 


Table  3.  Result  comparison  for  exercise  problems:  Heuristic  objective  value  (underlined,  if  different)  and  (W)-(H)  identical  objective  values. 


Goals 


Test 

case 

#  1  (Min): 
Unassigned 
task  parts 

#  4  (Min): 
Deviation 
level 

#  5  (Max): 
Spread 
primary 

#  6  (Max): 
Spread 
back-up 

#  7  (Min): 

M3  List 
position 

#  8  (Max): 
Residual 
salvo 

Solution  time 
(s):  Heuristic, 
(W),(H) 

R-Sl 

0 

6.5 

3 

2,4 

97 

148,  202 

1,  1,2 

R-S2 

0 

1.5 

4 

2,3 

132 

117,  184 

1,  1,2 

R-S3 

0 

49.0,  43.5 

8 

6 

481 

168,  257 

3,2,8 

R-S4 

0 

34.5 

6 

0 

151 

143,  245 

2,  1,4 

T-Sl-1 

0 

40.0 

4 

6 

330 

210,  253 

3,  3,9 

T-Sl-2 

0 

4.0 

4 

4,6 

540,  536 

174,212 

4,  10,  30 

T-Sl-3 

0 

1.7 

6 

4 

747,  721 

116,  166 

3,  11,47 

T-S2-4 

0 

46.7 

4 

6 

296 

233,  267 

4,  1,7 

T-S2-5 

0 

12.7 

4 

6 

502,  491 

198,  208 

7,  11,40 

T-S2-6 

0 

7.7,  2.7 

6 

6 

736,  684 

132,  179 

8,41,96 

T-S3-7 

0 

33.3 

4 

6 

364,  360 

203,  252 

3,2,  11 

T-S3-8 

0 

40.0 

6 

0,3 

584,  574 

151,  187 

6,  5,  39 

T-S3-9 

0 

2.7 

6 

0 

759,  718 

94,  123 

7,  29,  38 
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All  (W)  and  (H)  solutions  are  the  same  for  each  goal  in  all 
of  our  exercise  problems.  Note  that  goals  2  and  3  are  excluded 
because  exercise  problems  contain  neither  “penalty”  nor 
“expend”  firing  units.  As  opposed  to  the  lab  cases  where  goal 
4  was  commonly  the  highest  priority  objective  for  which  the 
heuristic  failed  to  produce  an  optimal  solution,  we  observe 
that  this  phenomenon  only  occurs  here  in  two  of  the  13 
cases.  For  the  remaining  cases,  it  is  either  goal  6,  7,  or  8 
which  returns  suboptimal  solutions.  The  overall  solutions 
(considering  all  objectives)  are  of  acceptable  quality.  Regard¬ 
ing  each  individual  goal,  for  most  cases,  the  heuristic  results 
are  reasonably  comparable  to  those  produced  by  the  exact 
models. 

However,  it  is  when  compared  to  the  “manual”  solution 
produced  by  the  TSC  (details  not  shown)  for  these  exercise 
cases  that  we  realize  the  benefit  of  using  optimization-based 
allocation.  Specifically,  the  TSC  fails  to  meet  goal  1  in  two  of 
the  four  Rimpac  cases  (leaving  as  many  as  eleven  unassigned 
task-parts  in  one  case),  and  in  eight  of  the  nine  Texas-Thunder 
scenarios  (leaving  up  to  four  unassigned  task-parts).  Com¬ 
pared  to  all  exact  solutions  (and  almost  all  heuristic  solutions) 
for  other  goals,  the  TSC’s  solutions  are  always  worse.  That 
is,  the  detriment  in  goal  1  (or  others)  is  not  incurred  for  the 
benefit  of  less  important  goals  either  in  those  solutions.  For 
example,  TSC’s  solution  gaps  with  respect  to  exact  solutions 
are,  on  the  average,  460%,  31%,  49%,  5%,  and  29%  for  goals 
4,  5,  6,  7,  and  8,  respectively. 

5.4.  Additional  Exercise-Based  Problems 

Three  other  examples  drawn  from  the  exercise  problems 
are  considered  as  separate  test  cases.  These  are  created  by 
including  all  salvos  within  the  same  Texas-Thunder  strike. 
That  is,  anew  “T-Sl”  case  combines  “T-Sl-1 ,”  “T-Sl-2,”  and 
“T-Sl-3”  salvos,  a  new  “T-S2”  case  comprises  “T-S2-4 ,”  “T- 
S2-5 ,”  and  “T-S2-6”  salvos,  and  a  new  “T-S3”  case  combines 
“T-S3-7 ,”  “T-S3-8 ,”  and  “T-S3-9.”  These  cases  have,  respec¬ 
tively,  168,  156,  and  175  tasks,  and  232,  229,  and  223  task- 
parts.  (That  is,  these  are  fleet- wide,  complex,  very  large-scale 
launch  plans.)  Because  all  task-parts  are  allocated  simultane¬ 
ously  for  all  salvos  in  a  strike,  the  “best-case  strike  planning” 
assumption  of  the  exercise  cases  is  dropped.  That  is,  a  mis¬ 
sile  allocated  to  a  redundant  task-part  in  one  salvo  cannot  be 
allocated  to  another  task-part  for  any  salvo  in  the  same  strike. 
This  is  conservative  under  the  assumption  that  redundant 
allocations  may  need  to  be  used  if  primary  missiles  fail. 

(W)  and  (H)  runs  for  “T-Sl ,”  “T-S2,”  and  “T-S3”  take 
considerably  longer  than  those  when  salvos  are  individually 
considered,  requiring  several  hours  of  computational  time. 
While  in  some  cases  the  TSC  may  have  sufficient  time  to 
undertake  such  analysis,  we  conservatively  set  up  these  cases 
to  run  in  a  maximum  of  10  minutes  or  to  a  relative  inte¬ 
grality  tolerance  of  1%,  whichever  comes  first.  For  the  (W) 


formulation,  this  set-up  needs  no  explanation.  For  the  (H^) 
formulation,  we  allow  up  to  5  minutes  for  k  =  1.  We  then 
allow  up  to  50%  of  the  remaining  time  for  k  =  2,  and  so  forth. 
The  (H,0  formulation  optimizes  (or,  at  least,  gives  a  bound) 
on  each  individual  objective  function  after  restrictions  are 
imposed  by  the  higher-priority  objectives.  This  allows  us  to 
determine  not  only  the  quality  of  the  (H^)  solution  but  also 
that  of  the  (W)  solution  for  that  goal. 

Both  exact  methods  typically  exhaust  the  allotted  time 
although  there  are  exceptions.  For  example,  the  (W) 
problem  may  terminate  early  within  1%  from  optimal  of  its 
own  weighted  sum  objective,  which  does  not  imply  1%  opti¬ 
mal  for  each  goal  individually  considered.  On  the  other  hand, 
the  (H)  solution  may  exhaust  the  maximum  time  for  one  of 
the  models  (H^),  but  terminate  early  overall  if  subsequent 
models  (H^),  for  k'  >  k ,  finish  within  1%  of  tolerance.  The 
Heuristic  algorithm  requires  between  15  and  42  seconds  for 
these  problems. 

In  all  three  new  cases,  solutions  for  goal  1  are  optimal,  leav¬ 
ing  4,  3,  and  21  nonallocated  task-parts,  respectively,  using 
either  formulation.  Notice  that  this  outcome  for  goal  1  con¬ 
trasts  with  that  of  all  other  cases  in  which  all  taskings  are  met. 
For  goal  4,  the  (H4)  model  returns  a  feasible  solution  in  each 
case,  but  the  bound  is  always  much  lower  than  such  solu¬ 
tion’s  objective  value.  For  example,  for  “T-Sl”  the  achieved 
value  is  1 .7  with  a  lower  bound  of  0.0.  Given  the  goal’s  nature 
and  the  achieved  value,  the  relative  error  is  of  little  impor¬ 
tance  in  practice:  after  the  allocation,  missile  inventories  left 
across  firing  units  are  very  well  balanced.  The  (W)  solution 
for  goal  4  is  always  the  same  as  or  worse  than  that  of  the  (H4) 
model.  For  all  remaining  goals  (5-8),  optimal  solutions  are 
achieved  by  the  (H^)  models  (with  similar  values  found  by 
the  (W)  model),  but  of  course,  those  solutions  could  change 
if  a  better  solution  for  goal  4  were  discovered.  The  heuris¬ 
tic  performance  is,  again,  of  high  quality,  yielding  optimal 
solutions  for  goal  1  in  two  of  the  three  cases. 

6.  IMPLEMENTATION 

To  demonstrate  our  models  and  methods,  and  to  gain  the 
trust  and  acceptance  of  navy  leadership,  we  had  to  build 
a  disposable  prototypic  graphical  user  interface.  However, 
we  knew  that  if  we  succeeded  in  gaining  acceptance  of 
our  model,  our  methods  would  be  integrated  into  a  combat 
system,  and  our  interface  would  be  discarded. 

Figure  4  shows  one  of  the  planning  screens  of  our  proto¬ 
type  system  used  in  tests  and  initial  implementation:  the  user 
interface  view  of  Tomahawk  task  allocations  to  firing  units. 
At  the  upper  right  are  navigation  buttons  to  other  screens. 
The  next  row  of  buttons  permits  scenario  import,  export, 
display  filtering,  printing,  and  emailing  of  a  scenario,  while 
the  clear  button  erases  a  candidate  solution,  and  the  allo¬ 
cate  button  runs  the  optimization.  For  each  task,  a  row  shows 
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Figure  4.  User  Interface  View  of  Tomahawk  Task  Allocations  to  Firing  Units. 


the  feasible  launch  area(s),  target  number,  mission  id,  task, 
day,  and  time.  The  required  number  of  primary,  ready- spare, 
and  backup  missiles  appear  in  columns  labeled  P,  RS,  and 
BU,  respectively.  The  column  entry  “Req”  specifies  the  pre¬ 
ferred  Tomahawk  variant  missile  to  be  used.  The  next  three, 
highlighted,  columns  show  the  allocation  of  the  required  task- 
parts  to  firing  units.  The  TSC  can  optimize  a  fleet-wide  or 
theater-wide  Tomahawk  strike  using  a  single  button  on  this 
dashboard  with  a  laptop  computer.  (Remark:  At  the  time 
the  prototype  was  being  developed  and  shown  to  the  navy, 
refinements  to  the  exact  formulation  for  better  solvability 
were  still  underway,  so  only  the  heuristic  algorithm  was 
implemented  within  the  graphical  user  interface.) 

We  used  this  interface  to  demonstrate  our  planning  model 
to  a  number  of  senior  officers,  and  in  presentations  to  larger 
groups,  such  as  in  meetings  of  the  Military  Operations 
Research  Society.  Subsequently,  the  U.S.  Navy  hired  a  con¬ 
tractor  who  efficiently  implemented  our  planner  in  the  opera¬ 
tional  environment,  and  thus  our  prototype  user  interface  was 
discarded. 

7.  CONCLUSION 

A  significant  part  of  our  research  has  been  devoted  to  sys¬ 
tem  analyses:  capturing  all  the  engineering  details  governing 
the  capabilities  and  constraints  associated  with  the  allocation 
of  Tomahawk  tasks  to  firing  units,  merging  the  engineer¬ 
ing  with  naval  tactics,  and  isolating  key  issues.  We  have 
also  developed  and  defined  a  standard  lexicon  with  which 
to  discuss  these  results. 

Historically,  a  TSC  primarily  performed  the  allocation  of 
tasks  to  firing  units  by  hand.  Our  challenge  was  to  replace 
the  manual  procedures  with  a  reliable,  carefully  documented 
system  to  rapidly  produce  high-quality  plans.  We  have  shown 
by  example  how  to  plan  optimally  for  a  variety  of  firing  units 
acting  in  unison  as  a  combat  strike  force.  We  now  know 
how  to  deliver  multiple- ship,  multiple-missile  strike  plans, 
taking  into  account  a  complex  variety  of  constraints  and 
objectives,  within  seconds. 

APPENDIX 

This  appendix  describes  our  simplified  model  formulation  for  AEGIS 
ships  only. 


Indices 


s  e  S  set  of  AEGIS  ship  firing  units  that  can  be  allocated 
task-parts 

h  e  H  set  of  VLS  half-modules  on  the  AEGIS  ship  firing  units 
that  can  be  allocated  task-parts 

t  e  T  set  of  strike  package  tasks  that  are  to  be  allocated  to  AEGIS 
ship  firing  units 

p  e  P  set  of  task-parts  (primary,  ready-spare,  backup) 
w  e  W  set  of  Tomahawk  weapon  types  (i.e.,  XMIDs)  loaded 
on  the  AEGIS  ship  firing  units  that  can  be  allocated  to 
task-parts 

i  G  I  set  of  Tomahawk  preparation  and  launch  time  interval 
conflict  task  sets  (constructed  based  on  the  preparation 
and  launch  time  intervals  of  the  Tomahawk  MSN,  and 
the  AEGIS  ship  firing  unit  restriction  that  only  one  Tom¬ 
ahawk  at  a  time  can  be  powered  up  in  a  VLS  half 
module) 

k  e  K  set  of  objective  functions,  {1, . . . ,  8} 
m  e  M,p  ordinal  set  for  the  Tomahawk  missiles  required  for 
task-part  p  of  task  t  (e.g.,  Mtp  =  0  if  no  missiles  are 
required,  or  Mtp  =  {1,2}  if  two  missiles  are  required  for 
the  task-part) 


Index  Maps  and  Subsets 


Lw  —  T 
Lwi  —  TSXJ0 

pNS  g  p 

Wsh  c  w 


wsht  c  Wsh 


Hs  c  H 


St£S 
SE  OS 
SFR  c  S 
SBU  c  S 


cTSC  cTSC 
A  ’  ° tp 


c  S 


subset  of  tasks  that  are  geofeasible  for  attack  by  ship 
3  with  Tomahawk  weapon  type  w 
subset  of  tasks  in  conflict  task  set  i  that  are  geofeasi¬ 
ble  for  attack  by  ship  3  with  Tomahawk  weapon  type 
w,  i.e.,  Tswi  —  Tsw  n  i 
subset  of  nonseparable  task-parts  for  task  t 
subset  of  Tomahawk  weapon  types  that  are  on  ship  3 
in  half-module  h 

subset  of  Tomahawk  weapon  types  that  are  on  ship  3 
in  half-module  h  that  can  be  used  for  task  t 
subset  of  AEGIS  ship  VLS  half-modules  that  are  on 
ship  5 

subset  of  AEGIS  ship  firing  units  that  are  geofeasibly 
located  to  perform  task  t 

subset  of  AEGIS  ship  firing  units  that  have  been 
designated  as  “expend”  ships 
subset  of  AEGIS  ship  firing  units  to  which  we  wish 
to  allocate  primary  task-parts  (see  Objective  5) 
subset  of  AEGIS  ship  firing  units  to  which  we  wish 
to  allocate  backup  task-parts  (see  Objective  6) 
subsets  of  AEGIS  ship  firing  units  designated  by 
the  TSC  as  candidates  to  which  to  allocate  task  t 
or  task-part  p  in  task  t,  respectively. 


Naval  Research  Logistics  DOI  10.1002/nav 


Newman  et  al.:  Optimizing  Assignment  of  Tomahawk  Cruise  Missiles 


293 


Data 

These  parameters  are  needed  only  for  the  objective  function  definitions. 
ps  penalty  for  using  ship  s 

m]wp  penalty  for  allocating  weapon  w  to  part  p  in  task  t  (based  on  M3 
position  of  weapon  w  for  task  t,  and  the  weight  of  part  p  to  match 
the  M3  list) 

vw  relative  value  of  weapon  w  for  residual  salvo  computation 
nshw  number  of  weapons  of  type  w  in  half-module  h  on  ship  s 
cok  weight  given  to  goal  k  when  a  weighed  objective  is  considered 

Primary  Decision  Variables 

Xshwtpm  equals  1  if  ship  s  is  assigned  to  perform  the  ra-th  missile 
requirement  for  task-part  p  in  task  t  (tentatively  using  a 
weapon  of  type  w  in  half-module  h ),  and  0  otherwise 
Ystpm  equals  1  if  ship  s  is  assigned  to  perform  the  ra-th  missile 
requirement  for  task-part  p  in  task  t,  and  0  otherwise 
Ufpm  equals  1  if  the  ra  -th  missile  requirement  for  task-part  p  in  task 

t  is  unmet,  and  0  otherwise 

Auxiliary  Decision  Variables 

These  variables  are  used  only  to  define  the  objective  functions. 

Zk  value  of  objective  function  k 
Vs  equals  1  if  ship  5  is  used,  and  0  otherwise 
RD  S  difference  between  the  residual  number  of  missiles  on  ship 

^  and  the  average  residual  number  of  missiles  across  all 
“non-expend”  ships 

PR  V  equals  1  if  any  primary  task-part  is  allocated  to  ship  s  (only 
for  5  e  SFR),  and  0  otherwise 

BU5  equals  1  if  any  back-up  task-part  is  allocated  to  ship  s  (only 
for  ^  £  SBU),  and  0  otherwise 

NS stp  equals  1  if  all  missile  requirements  for  nonseparable 

task-part  p  of  task  t  are  assigned  to  ship  5,  and  0  otherwise 
Dshw  equals  1  if,  based  on  the  information  held  by  the  TSC,  one 
weapon  of  type  w  would  be  the  most  valuable  weapon  to 
remain  on  ship  s  in  half-module  h  after  firing  all  primary 
task-part  missiles,  and  0  otherwise 


Objective  Functions 


Objective  k  =  1: 

min  z\  — 

Y  Ufpm 

teT,peP, 

meMtp 

Objective  k  =  2: 

min  Z2  — 

Ep°v* 

seS 

Objective  k  =  3: 

max  Z3  — 

^  '  Ystpm 

t  eT  ,s  eStC\SE , 
peP,meMtp 

Objective  k  =  4: 

min  Z4  — 

E  RDv 

seS\SE 

Objective  k  =  5: 

max  Z5  = 

E  PR> 

seSFR 

Objective  k  =  6: 

max  z.6  — 

E  BU' 

Objective  k  =  1: 

min  zi  — 

m?twpXshwtpm 

t  eT  ,s  eSt  ,heHs , 
weWsht,peP, meMtp 

Objective  k  =  8: 

max  zg  = 

E  vwDshw 

seS,heHs,weWsh 


Constraints  Enforcing  Physical,  Policy,  and  whim  of 
TSC  Requirements 

Ensure  weapon  availability  for  tasks  allocated  to  half-modules: 

^  '  Xshwtpm  —  ^shw  Vs  £  S,h  £  HS,W  £  Wsh  (1) 

t^Tsw,peP  ,meMtp 

Establish  whether  or  not  a  specific  firing  unit  receives  a  task  allocation: 

Xshwtpm  —  Ystpm  Vt  £T,S£St,p£P,m£  Mtp  (2) 

heHs,weWsht 

Preclude  multiple  assignments  of  the  same  task: 

E  YsfPm  -  1  Vt  £  T ,  p  £  P ,m  £  Mtp  (3) 

seSf 

Preclude  partial  completions  of  task  parts: 

Y  Yst,' primary7/ 1'  =  YstPm  Vt  £  T ,  p  £  P ,m  £  Mtp  (4) 
seSt  seSt 

Do  not  assign  more  than  one  conflicting  task  to  the  same  half-module: 

Xshwtpm  <  1  Vs  £  S,h  £  Hs,i  £  I  (5) 

weWsh,teTSWi, 
peP  ,meMtp 

Allocate  each  task,  or  signal  a  violation: 

Ystpm  +  Utpm  —  1  Vt  £  T ,  p  £  P ,m  £  Mtp  (6) 

seSf 

Enforce  primary  and  back-up  task-part  relations: 

V  _j_  V  <<  j 

1  st, 'primary  ,m  '  1  sf/back— up' ,m'  —  1 

Vt  £  T,St  £  S,m  £  M primary',  m'  £  Af^/back_up/  (7) 
Enforce  identical  ship  allocation  for  nonseparable  task-parts: 

E  Y°>p"  =  \M<p\  NS*>p  Vf  e  T’s  e  s»p e  P>N  w 

meMtp 

Enforce  primary  and  ready- spare  task-part  relations: 

^st, 'primary', m  t,' ready— spare' ,m  ^  ^  T,S  £  Sf,tTl  £  47^/j-eady—spare'  (9) 

Exclude  any  allocation  the  TSC  dislikes: 

Ystpm  =  0  Vt  e  T,  p  e  P,m  e  MIp,s  £  sjsc  n  sjpc  (10) 
Binary  decision  variables: 

Xshwtpm  ^  {0,  1}  Vt  £  T ,S  £  St,h  £  HS,W  £  Wshtim  £  Mtp  (11) 

NSstp  €  {0, 1)  Vt  eT,s  e  S,,p  e  P,N  (12) 

Continuous  decision  variables: 

Ystpm  e  [0, 1]  Vt  £T,S  £  st,p  £  P,m  £  Mp  (13) 

Utpm  0,1]  Vt  £  T ,  p  £  P ,m  £  Mtp  (14) 
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Constraints  Defining  Objective  Functions 

These  relationships  do  not  explicitly  constrain  the  solution  and  may  be 
handled  implicitly  in  some  implementations.  They  can,  of  course,  be  omitted 
if  the  corresponding  objective  is  not  considered  in  a  particular  scenario. 

For  objective  k  =  2: 

Ystpm  —  14,  Vt£T,sESt,p€zP,m€:Mp 

V5e{0,l},  Vs  g  S 

For  objective  k  =  4: 

RD,  >0,  Vs  e  S 

There  are  additional  details  given  in  Demir  [9]. 

For  objectives  k  =  5  and  k  =  6: 

PRS  <  J2  Vprimaiy',m>  €  SPR 
Vt\seSt, 

mGMt /primary' 

BUj  <  Y.  Y«! back— up',m-  Vi  6  SBU 

Vt\seSt, 
mGMt'  back-up' 

0  <  PR,  <1,  Vs  e  SPR 
0  <  BU,  <  1,  Vs  €  SBU 

For  objective  k  =  8: 

^shw  ^  '  X shwt, 'primarj/ ,m  Dshwi  Y?  £  S,h  G  Hs,  W  G  Wsh 

teTsw,m£  Mt /primary/ 

£  Dshw<^  Vies,  heHs 

WGWsh 

0  <  Dshw  <1,  Vs  G  S,  h  G  Hs,  W  G  Wjfc 
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